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ABSTRACT 

We present evidence for a strong correlation between the concentration of bulges and the mass of 
their central supermassive black hole (Mbh) — more concentrated bulges have more massive black holes. 
Using C re (l/3) from Trujillo, Graham, & Caon (2001b) as a measure of bulge concentration, we find that 
log(M b h/M ) = 6.81 (±0.95)C re (1/3) + 5.03 ±0.41. This correlation is shown to be marginally stronger 
(Spearman's r s = 0.91) than the relationship between the logarithm of the stellar velocity dispersion 
and logMbh (Spearman's r s — 0.86), and has comparable, or less, scatter (0.31 dex in logMbh, which 
decreases to 0.19 dex when we use only those galaxies whose supermassive black hole's radius of influence 
is resolved and remove one well understood outlying data point). It would appear that the central black 
hole mass can be estimated from surface photometry alone, without the expensive addition of velocity 
dispersion determinations. 

Subject headings: black hole physics, galaxies: fundamental parameters, galaxies: kinematics and 
dynamics, galaxies: nuclei, galaxies: photometry, galaxies: structure 



1. INTRODUCTION 

Observations over the last few years have established 
that supermassive black holes (SMBHs; ~ 10 6 - 1O 9 M ) 
are a common, if not universal, feature at the centers of 
elliptical galaxies and the bulges of early-type spirals (Ko- 
rmendy & Richstone 1995; Magorrian et al. 1998). Initial 
correlations between the masses of SMBHs and the abso- 
lute B-band luminosities of the host bulges 3 had a large 
scatter (~ 0.5 — 0.6 dex in logAfbh, but see McLure & 
Dunlop 2001) which could not be accounted for by the as- 
sumed measurement errors. Subsequent studies uncovered 
a stronger correlation between the mass of the SMBH and 
the stellar velocity dispersion of the bulge, with consider- 
ably smaller scatter: only ~0.3 dex in log Mbh (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). Merritt & Ferrarese 
(2001a) argued that the observed scatter, for a select sam- 
ple of 12 galaxies thought to have the most reliable SMBH 
mass estimates, was fully consistent with the measurement 
errors alone: that is, there may be no intrinsic scatter in 
the correlation. This clearly suggests that a strong cross- 
talk exists — or once existed — between the central black 
hole and its host bulge. The reasons for this, and the pre- 
sumed formation mechanism are, however, not well under- 
stood, although many possibilities have been put forward 
(see, e.g., the discussion in Merritt & Ferrarese 2001b). 

Recently, Graham, Trujillo, & Caon (2001) have shown 
that the central concentration of bulge light, measured 
within one effective half-light radius, positively correlates 
with the stellar velocity dispersion of the bulge. Following 
up on this, we present here, for the first time, a correlation 
between SMBH mass and bulge concentration, showing 
that more concentrated bulges have more massive SMBHs. 
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3 By the term bulge we mean both elliptical galaxies and the bulges 

4 For three galaxies, we used HST F814W images instead. 



This correlation is found to be at least as strong as that 
between SMBH mass and stellar velocity dispersion, and 
may have less scatter. We suggest that the bulge con- 
centration is certainly as significant a parameter, and one 
perhaps more revealing, than velocity dispersion (which 
presumably is a response to the underlying bulge mass 
distribution) for understanding the fueling and growth of 
central SMBHs. Furthermore, bulge concentration is eas- 
ier to measure. 



2. GALAXY DATA AND MEASUREMENTS 

We began with the updated list of SMBH detections and 
mass estimates in the first two sections of Merritt & Fer- 
rarese's (2001b) Table 1. These black hole masses include 
a number of revised estimates from the "Nuker group" 
and STIS IDT team presented by Kormendy & Gebhardt 
(2001). This is an initial total of 30 galaxies, including the 
Milky Way. The only quantity that we have changed is the 
SMBH mass estimate for NGC 4374. Although this object 
appears in the list of galaxies with "reliable" SMBH mass 
estimates (because the black hole's sphere of influence has 
been resolved), Maciejewski & Binney (2001) recently re- 
vised its mass estimate lower by a factor of four, after 
taking proper account of finite slit-width effects. 

We searched the various public archives for high-quality 
i?-band images 4 which were large enough to guarantee 
good sky subtraction and which had no central saturation. 
We found useful images for a total of 21 galaxies, primar- 
ily from the Isaac Newton Group and Hubble Space Tele- 
scope (HST) Archives; we also used images from Frei et al. 
(1996) and Hintzen et al. (1993), available via the NASA 
Extragalactic Database (NED). We were also able to use 



of spiral galaxies. 
1 



2 



GRAHAM, ERWIN, CAON, TRUJILLO 



an unpublished image obtained with the WIYN Telescope 5 
for NGC 3245, and the minor-axis near-infrared surface- 
brightness profile of the Milky Way published by Kent, 
Dame, & Fazio (1991), making a total of 23 galaxies with 
useable data. 

A full discussion of the images for individual galaxies, 
along with reduction procedures and the extracted light 
profiles, analysis and model fitting for each galaxy will 
be presented in Erwin et al. (2001). Briefly, we fitted el- 
lipses to the isophotcs, allowing the position angle and 
cllipticity to vary with radius. The resulting light pro- 
files were then fitted with a seeing-convolved 6 Sersic (1968) 
j- 1 /™ model. We modeled disk galaxy profiles with a com- 
bined (seeing-convolved) exponential disk and r 1 '™ bulge. 
For two galaxies with strong bars, we used the light profile 
derived from cuts along an axis perpendicular to the bar; 
this produced much better fits and avoided most of the 
bar contribution. The inner arc second of the profiles was 
generally excluded from the fit, since these regions are of- 
ten dominated by relatively flat cores, bright nuclear disks, 
or nuclear point-sources (e.g., Rest et al. 2001, and refer- 
ences therein), none of which can be modeled with Sersic 
profiles. We are thus measuring the overall concentration 
of the bulge, independent of any separate central stellar 
components like nuclear disks. Only two galaxies could 
not be successfully modeled. The final set of 21 galaxies 
with well- fitted profiles is given in Table 1. 

We then computed the concentration of the best-fitting 
r 1 /™ models using the central concentration index first pre- 
sented in Trujillo, Graham & Caon (2001b) and further 
developed in Graham et al. (2001). This index measures 
the light concentration within a bulge's half-light radius 
(r e ): it is the ratio of flux inside some fraction a of the 
half-light radius to the total flux inside the half-light ra- 
dius. For an r 1 /™ model, this index can be analytically 
defined as 

where n is the shape parameter of the r 1 /™ model and b n is 
derived numerically from the expression r(2n)=27(2n, b n ) 
where T(a) and "f{a,x) are respectively the gamma func- 
tion and the incomplete gamma function (Abramowitz & 
Stegun 1964). (This index can be measured empirically, 
without the use of Sersic models, but for bulges in disk 
galaxies this would first require successful two-dimensional 
modeling and subtraction of disks, bars, etc.) The pa- 
rameter a can be any value between and 1, and defines 
what level of concentration is being measured. Follow- 
ing Graham et al. (2001), we used a value of a = 1/3. 
We did however explore a range of values, finding that 
a = 1/3 roughly produced the minimum vertical scatter in 
the logMbh — C re (a) correlation. CV c (l/3) is then simply 
the ratio of flux inside one-third of the half-light radius 
to the flux inside the entire half-light radius, (which is 
of course half the total bulge luminosity). The CV e (l/3) 

5 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the National 
Optical Astronomy Observatories. 

6 We used a Moffat function to describe the point-spread function; seeing full-width half-maxima were measured from stars in the individual 
images. 

7 This issue is discussed at length in Merritt & Ferrarese (2001a) for one galaxy in particular, namely, the Milky Way. 

8 These data points may be revealing negative curvature in the log M^-C re (1/3) relation, but wc consider it premature to reach such a 
conclusion based on one data point at each end of the relation. 



values are listed in the final column of Table 1. Because 
these values are analytically derived from the best-fitting 
Sersic index n, the uncertainty in CV e (l/3) depends di- 
rectly on the uncertainty in n and is derived by standard 
propagation of errors. Error estimates for n are based on 
the results of Caon, Capaccioli, & D'Onofrio (1993), who 
found a typical uncertainty of ~ 25% when fitting with 
Sersic profiles. For Sersic values of n between 2 and 11, 
this corresponds to a 10-15% uncertainty in the bulge con- 
centration. 

For comparison with the known SMBH mass-velocity 
dispersion relation, we also list the velocity dispersions a c 
and corresponding errors for each galaxy; these are taken 
from Merritt & Ferrarese (2001b) and thus incorporate 
the equivalent-aperture correction described in Ferrarese 
& Merritt (2000). As Merritt & Ferrarese (2001a) showed, 
these values differ, on average, by only ~ 1% from the a e 
values used by Gebhardt et al. (2000). 

3. RESULTS 

Correlations between SMBH mass and bulge concentra- 
tion are presented in Figure 1; for comparison, we also 
show the correlations between SMBH mass and velocity 
dispersion for the same galaxies. We used the bisector 
linear-regression routine from Akritas & Bershady (1996) 
to fit a line to these correlations. This regression routine 
treats both variables equally, and allows for intrinsic scat- 
ter as well as measurement errors in the data; as Merritt & 
Ferrarese (2001a) point out, it is generally the best method 
to use when there are errors in both variables. Using 
the orthogonal regression analysis of Akritas & Bershady 
(1996) and the orthogonal distance regression routine FI- 
TEXY of Press et al. (1989) — alternate methods for data 
sets with errors in both variables — gave consistent re- 
sults. We computed the Pearson correlation coefficient r 
and Spearman rank-order correlation coefficient r s , both 
of which are given in Figure 1 . The Spearman coefficient is 
preferred as it is more robust to outliers and does not pre- 
suppose a linear relation. The best linear fit to the whole 
sample is logM bh = 6.81(±0.95)CV e (l/3) + 5.03 ± 0.41. 

Figure 1 shows that the correlation between black hole 
mass and bulge concentration is extremely good — as good 
as or better than that between black hole mass and ve- 
locity dispersion. In addition, the low x 2 value of 0.82 
suggests a scatter consistent with the measurement errors 
alone, implying negligible intrinsic scatter (as Ferrarese & 
Merritt 2000 claimed for the SMBH - velocity dispersion 
relation). This conclusion does, however, depend on how 
well determined the errors are; see Erwin et al. 2001. 

Data points at the extreme ends of a correlation can 
be very useful for determining the true slope, due to the 
weight they lend, but by the same token they can heavily 
bias the data to produce a misleading slope if they them- 
selves have not been well determined 7 . We have identified 
two such potential outliers 8 (NGC 7457 and NGC 6251) in 
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Figure la (see Section 4). We have repeated our analysis 
without them (bottom panels of Figure 1). In so doing, the 
log Mbh-log a c relation is even weaker than the logMbh- 
C rc (1/3) relation; it also has 27% more scatter in logMbh- 
The list of galaxies in Merritt & Ferrarese (2001b), from 
which we constructed our sample, was subdivided accord- 
ing to whether or not the black hole's sphere of influence 
had been resolved 9 Of the twenty-two "resolved" galaxies, 
we have images and useful fits for fourteen. We rederived 
the relations seen in Figure 1 using this smaller sample, 
and found that the strength of both correlations improved; 
for this subsample, both correlations appear equally strong 
(Figure 2). 

A word of caution may be in order when comparing dif- 
ferent measures of significance for these relations. The 
strength of a correlation itself — regardless of which func- 
tion fits it — is best measured by the Spearman rank- 
order coefficient r s . The \ 2 merit function for a linear 
fit to the data, the Pearson coefficient r, and the verti- 
cal scatter in log Mbh all measure how well a straight line 
fits the data (or the logarithm of the data, as the case 
may be). In this vein, we note that while the strength of 
the log Mbh _ CV e (a) correlation — as measured with the 
Spearman rank-order correlation coefficient — remains un- 
changed as a and hence C rc (a) vary, the x 2 merit function 
steadily, and significantly, decreases as a increases. This 
is easily understood from the way the absolute error in n 
propagates to an absolute error in C re (a). Furthermore, 

The x 2 value depends on the size of the measurement er- 
rors: overestimating the errors will decrease the resulting 
X 2 , even though the correlation is unchanged; underesti- 
mating the errors can produce a mislcadingly large x 2 ■ 
Thus, even though the x 2 values for the M^-C Tc (1/3) re- 
lation are either the same as (Figure 2) or much smaller 
than (Figure 1) those for the M^-a c relation, we do not 
take that as strong evidence that the Mbh-C rc (l/3) rela- 
tion is better. 

Ferrarese & Merritt (2000) argued that their optimal 
logMbh — logcr c relation had negligible intrinsic scatter; 
this led them to posit that, "Our results suggest that the 
stellar velocity dispersion may be the fundamental param- 
eter regulating the evolution of supermassive BHs in galax- 
ies." All twelve of their 'Sample A' galaxies, from which 
this conclusion was reached, had an uncertainty of ±13% 
on their central velocity dispersions, except for the Milky 
Way which had an uncertainty of ±20%. If these uncer- 
tainties have been over-estimated it will result in an un- 
derestimate to the x 2 value of the fit which may then lead 
one to wrongly conclude that there is no intrinsic scatter 
in the relation. The situation is identical if the ~10-15% 
errors we assigned to the central concentration indices are 
to large. 

To investigate the effect that decreasing the error esti- 
mates has on our correlations, we reperformed the regres- 
sion analysis assuming only a 10% error on n, which trans- 
lates to an impressively small 3% error in the derived con- 
centration index when n — 8 and a 6% error when n = 1. 
For the full galaxy sample, the slope of the best-fitting 



line in the logMbh _ CV e (1/3) diagram decreased slightly to 
6.12±0.78 and A log Mbh changed to 0.30 dex. Removing 
NGC 6251 and NGC 7457, the slope was 6.04±0.53 and 
A log Mbh decreased to 0.24 dex. 

The bulges studied here clearly have different luminos- 
ity distributions and, unless M/L varies with radius in a 
very contrived fashion, they also have different mass dis- 
tributions from each other. If they did all possess the same 
universal structure, then the concentration index would be 
constant and identical for every bulge, which it is not. We 
will never have an accurate picture of bulge formation if 
we continue to pigeon-hole bulges into two simple classes, 
namely, r 1 / 4 and exponential. Graham (1998) wrote 
"...one might expect [Sersic's] n to correlate with the size of 
the black hole predicted to be at the center of many ellipti- 
cal galaxies." Since C r %l/3), as defined in equation 1, is a 
monotonic function of the global shape parameter n (Tru- 
jillo et al. 2001b), the SMBH mass - C re (l/3) correlation 
implies a correlation between SMBH mass and n as well. 
For the full galaxy sample, performing the bisector regres- 
sion analysis on logMbh and logn (assuming a 25% error 
in n) gives logM bh = 2.93(±0.43) logn + 6.42 ± 0.22 with 
a scatter of 0.32 dex in log M bh - Excluding NGC 7457 and 
NGC 6251 gives log M bh = 3.00(±0.17) logn + 6.45 ±0.11 
with a scatter of only 0.25 dex. The strength of these cor- 
relations are of course equal to those shown in Figure la 
and c. 

4. DISCUSSION 

Due to the weight that data points at the ends of a cor- 
relation can have on a line-of-best-fit, we have identified 
two outliers (NGC 7457 and NGC 6251) in Figure 1 which 
may be biasing the relation defined by the remaining bulk 
of data points. The most significant outlier in our relation 
is probably NGC 6251 in the top right of Figures la and 
2a; it has both the highest black hole mass and the high- 
est concentration index. There is reason to believe that its 
concentration index may be significantly in error. The best 
fitting Sersic r 1 /™ model to this galaxy has n = 11, which 
means the outer profile of this model declines slowly with 
radius; the observed light-profile only extends to 1 model 
half-light radius. The larger the value of n, the closer the 
r 1 /™ model approaches a power-law in behavior — having 
infinite brightness and an infinite half-light radius (e.g., 
Graham et al. 1996), resulting in excessively high values 
of C Tc (1/3). Values of n greater than about 10 should thus 
be treated with caution. The outlying point (NGC 7457) 
in the lower left of Figure la is less easily dismissed, and 
may be a true outlier worthy of individual investigation 10 
There is some evidence for a weak bar or lens in this galaxy 
(Michard & Marchal 1994). While we derived a light pro- 
file perpendicular to the major-axis of this feature, it does 
still have a finite width which can bias our Sersic fit to the 
bulge, giving a spuriously large value for n and hence an 
over-estimation of the bulge concentration. 

Although we cannot say which parameter (CV e (l/3) or 
logu) is better, we can identify some of the advantages 

9 This is not necessarily a guarantee of accurate mass estimates: NGC 4374, which is near the top of Merritt & Ferrarese's reliability ranking, 
recently had its mass readjusted by a factor of four (Maciejewski & Binney 2001). 

10 Interestingly, Tonry & Davis (1981) derived a stellar velocity dispersion for NGC 7457 of 129 km s _1 and Dressier & Sandage (1983) obtained 
a value of 136 km s" 1 . If these values are a better measurement than the value of 73 km s _1 which has been used, one can see that this galaxy 
would also be an outlier in Figure lb. 
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and disadvantages of each. Use of the concentration index 
for studies such as this may not be applicable to morpho- 
logically disturbed galaxies which may have undergone re- 
cent mergers or interactions (this could also influence the 
stellar velocity dispersion). Dominant cD galaxies that 
have acquired large extended envelopes can also be dif- 
ficult to model and, depending on the extent of the ac- 
creted material, their concentration index may be heavily 
biased. The stellar velocity dispersion may, on the other 
hand, be a more stable quantity in such cases. Velocity 
dispersions have additionally been measured for numerous 
(mostly nearby) galaxies. 

One obvious practical advantage that the concentration 
index has over stellar velocity dispersions is that images 
are far less expensive to acquire (in terms of telescope 
time) than stellar velocity dispersions; especially those at 
one effective radii. This is particularly important for stud- 
ies of high-redshift galaxies. Second, use of the bulge con- 
centration circumvents concerns that the SMBH may be 
influencing the luminosity-weighted nuclear velocity dis- 
persion measurements (e.g. Wandel 2001). Similarly, while 
n and C Te are global parameters, velocity dispersion mea- 
surements are affected by: kinematical sub-structure at 
the centers of bulges, rotational velocity, seeing conditions, 
and aperture-size. It should also be noted, however, that 
the presence of bars, rings, and lenses within spiral galaxies 
can complicate the determination of the bulge concentra- 
tion. A fourth advantage that the concentration index has 
is that it can be measured from photometrically uncali- 
brated images, it does not rely on distance measurements, 
rcdshift dependent corrections, or assumed mass-to-light 
ratios. 

It appear that the formation mcchanism(s) behind 
bulges must simultaneously determine their total lumi- 
nosity, their eventual luminous structure (as measured by 
concentration index and Sersic n) and mass distribution, 
their stellar velocity dispersion, and the final mass of the 
central SMBH. Subsequent evolution, including the effects 
of mergers (once this process has neared completion), ev- 
idently preserves the above connections. To date, most 
models incorporating SMBHs have addressed their forma- 
tion from the standpoint of the older Mbh--M"buige rela- 
tion (Haiman & Locb 1998; Richstone et al. 1998; Silk & 



Rees 1998; Blandford 1999; Kauffmann & Haehnelt 2000; 
Archibald et al. 2001; Umemura 2001); a few recent pa- 
pers have offered possible explanations for the M^-a cor- 
relation (Haehnelt & Kauffmann 2000; Adams, Graff, & 
Richstone 2001; Burkert & Silk 2001). We argue that a 
more complete understanding will be achieved when the 
correlation between SMBH mass and bulge concentration 
is also explained. 

More luminous (massive) bulges have larger values of 
n (Trujillo et al. 2001b, and references therein), greater 
central concentrations, deeper gravitational potential wells 
with steeper central gradients (Ciotti 1991, Trujillo et al. 
2001a), and more massive SMBHs. One might expect 
these characteristics to result in bulges which are better 
able to fuel and build their central black holes. It is how- 
ever likely that the processes which shaped the galaxy and 
built the supermassive black hole operated in tandem. 

Velocity dispersion measurements may simply be a 
somewhat more expensive tracer of the fundamental un- 
derlying mass distribution which can be more cheaply (in 
terms of telescope time) measured from galaxy light pro- 
files. 

We wish to thank Matthew Bershady for making avail- 
able the linear regression code from Akritas & Bershady 
(1996). 

Based on archival data obtained with the Isaac Newton 
Group of Telescopes operated on behalf of the UK Particle 
Physics and Astronomy Research Council (PPARC) and 
the Nederlanse Organisatie voor Wetenschappclijk Onder- 
zoek (NWO) on the island of La Palma at the Spanish 
Observatorio del Roque de Los Muchachos of the Insti- 
tuto de Astroffsica de Canarias. 

Based on observations made with the NASA/ESA Hub- 
ble Space Telescope, obtained from the data archive at 
the Space Telescope Institute. STScI is operated by the 
association of Universities for Research in Astronomy, Inc. 
under the NASA contract NAS 5-26555. 

This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. 
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Table 1 
Galaxy Parameters 



Galaxy 


Revised Hubble cr c 




n 


C re (l/3) 




Type (km s _1 ) 


(1O 8 M ) 





Ellipticals 



NGC 


821 


E6 


196 


± 


26 


NGC 


3377 


E5-6 


131 


± 


17 


NGC 


3379 


El 


201 


± 


26 


NGC 


4261 


E2-3 


290 


± 


38 


NGC 


4291 


E 


269 


± 


35 


NGC 


4374 


El 


286 


± 


37 


NGC 
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a NGC 2778 is classified as an elliptical galaxy in the NASA Extragalactic Database 
(NED), but its light profile clearly indicates an SO galaxy, with both an exponential 
disk and a central bulge; this interpretation is supported by the kinematical study of 
Rix, Carollo, & Freeman (1999). 

b HST F814W image 
c Thuan-Gunn r image 

d 2.4-/im minor-axis profile from Kent, Dame, & Fazio (1991). 



Note. — Galaxy types are from NED. Velocity dispersions and black hole masses 
from compilation in Merritt & Ferrarese (2001b), except SMBH mass for NGC 4374 
(revised by Maciejewski & Binney 2001; updated errors from Kormcndy & Gcbhardt 
2001). The central concentration C*r e (l/3) of each bulge was measured from ij-band 
images, except where otherwise noted, with a 25% uncertainty assumed for n. 
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Fig. 1. — Correlations between supermassive black hole mass and a) bulge concentration and b) stellar velocity dispersion within r e /8. 
c) and d) show the correlation after removing the two outlying galaxies (NGC 7457 and NGC 6251; see Section 3 and 4). The Milky Way 
(MW) has been indicated. The Spearman rank-order correlation coefficient r s is given, as is the Pearson linear correlation coefficient r. The 
X 2 merit function for a linear fit and the absolute vertical scatter A log Mbh are also given. Elliptical galaxies are denoted by filled circles, 
lcnticulars and spirals by open circles. 




^(1/3) a c (km s" 1 ) 

Fig. 2. — Same as Figure 1, except using only those galaxies with resolved SMBH spheres of influence (top section of Table 1, Merritt & 
Ferrarese 2001b). The dashed line shows the correlation after removing NGC 6251, as done in Figure 1. The slope to the dashed line in panel 
a) is 6.49±0.78 and has a vertical scatter of 0.19 dex. The slope to the dashed line in panel b) is 4.41±0.66 with a vertical scatter of 0.20 dex. 



